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We present a field-theoretic study of the critical Casimir force of the Ising universality class in 
a d- dimensional L't' 1 x L slab geometry with a finite aspect ratio p = L/L\\ above, at, and below 
T c . The result of a perturbation approach at fixed dimension d = 3 is presented that describes the 
dependence on the aspect ratio in the range p > 1/4. Our analytic result for the Casimir force 
scaling function for p — 1/4 agrees well with recent Monte Carlo data for the three-dimensional 
Ising model in slab geometry with periodic boundary conditions above, at, and below Te- 



la the presence of fluctuations with long-range correla- 
tions, so-called Casimir forces occur in macroscopic con- 
fined systems. The existence of such forces due to long- 
range critical fluctuations have been predicted by Fisher 
and de Gennes [jj for fluid films. For such systems with 
isotropic interactions, the critical Casimir forces depend 
only on the boundary conditions (b.c.) and on the ge- 
ometry of the confining surfaces as well as on the uni- 
versality class of the critical point 0, Q . For anisotropic 
systems (e.g., magnetic systems with noncubic symme- 
try), critical Casimir forces also depend on nonuniversal 
anisotropy parameters 0, Q . 

Considerable theoretical effort has been devoted to the 
study of critical Casimir forces in isotropic film systems 
over the past two decades @, 0, Q • In the present paper 
we shall focus on the Ising universality class withperiodic 
b.c. for which detailed Monte Carlo (MC) data (9|,E(| are 
available. While progress has been achieved by means of 
the e = 4 — d expansion [6j, |8| above the bulk critical 
temperature T c no theoretical prediction is available as 
of yet for the region below T c . The most interesting fea- 
ture is the existence of a pronounced minimum of the 
finite-size scaling function of the Casimir force below T c 
which is characteristic also for other film systems with 
realistic b.c. 0, EH EH, EH- In this Letter we present 
the result of a renormalization-group calculation within 
the framework of the </? 4 theory at fixed dimension d = 3 
0, [Hj] that is in good agreement with the MC data 0, E(| 
including the minimum below T c and the Casimir ampli- 
tude at T c . 

All of the existing theoretical studies [1, @, [1] of the 
critical Casimir force in film systems have considered an 
oo 2 x L geometry. This geometry is, of course, an ideal- 
ization that is only approximately realized in experiments 
or computer simulations. In fact, the MC simulations for 
the Ising universality class with periodic b.c. 9, 10] have 
been carried out for periodic L 2 x L slabs with finite as- 
pect ratios p = L/Lu in the range 1/14 < p < 1/3. Most 
of the available data are for p = 1/6. This appears to 
be well justified as the dependence on p for p <C 1 is 
expected to be rather weak. In Ref. Q it was stated 
explicitly that the MC results for p — 1/4 can hardly be 
distinguished from those for smaller values of p. 

Our new approach to the problem takes advantage of 



the fact that an L 2 x L finite-slab geometry is conceptu- 
ally simpler than a oo 2 x L film geometry for two funda- 
mental reasons. First, there exists no film transition at 
finite p > 0, thus there is no necessity of dealing with the 
as yet unsolved problem of dimensional crossover between 
the 3-dimensional bulk transition and the 2-dimensional 
film transition. Second, for p > 0, the system has a dis- 
crete mode spectrum with only one single lowest mode, 
in contrast to the more difficult situation of a lowest- 
mode continuum in film geometry. This opens up the op- 
portunity of building upon the advances that have been 
achieved in the description of finite-size effects in systems 
that are finite in all directions EH, EE EH ■ It is not 
clear a priori, however, in what range of p such a theory is 
reliable since, ultimately, for sufficiently small p <C 1, the 
concept of separating a single lowest mode must break 
down. Therefore, as a crucial part of our theory, we first 
provide quantitative evidence for the expected range of 
applicability of our theory at finite p. 

We start from the standard 0(n) symmetric isotropic 
Landau-Ginzburg- Wilson Hamiltonian 



H = Jd d r[^ + \{V^ 
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for the n component vector field ip(r) in a d-dimensional 
Ljj -1 x L finite-slab geometry with periodic b.c. in all 
directions. The fundamental quantity from which the 
critical Casimir force can be derived is the singular part 
f s (t, L, Ly) of the free energy / per unit volume and per 
component, divided by fc^T. The expected asymptotic 
(large L, large Ln, small t = T — T c ) finite-size scaling 
form of f s for isotropic systems is [17| 



f s (t,L,L ll )=L- d F(x,p) 



(2) 



with the scaling variable x = t(L/^o) 1 ^ v where £o is the 
amplitude of the bulk correlation length above T c . 

To study the p dependence of the scaling function 
F(x,p) we first consider the large-n limit (at fixed won). 
As an exact result we find in three dimensions 



F(x,p) = (8tt)- 
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with K(z) — X)m=-oo exp(— zm 2 ) where P(x,p) is de- 
termined implicitly by P = x — 4ir Qi (P 2 , p) . 

The amplitude F(0, p) at T c for 71 = oo in three dimen- 
sions is shown as thin solid line in Fig. 1. It interpolates 
smoothly between the limits of p = (film) and p = 1 
(cube). It is a monotonically decreasing function of p 
since the value of F (0, p) is suppressed as the confinement 
becomes stronger. As a nontrivial feature, Fig. 1 exhibits 
a negligible dependence on p for small p up to p < 1/4. 
The weak dependence of F(x,p) for p < 1/4 also per- 
tains to the central finite-size region |x| < O(l) around 
T c . This suggests that studying F(x,p) in a finite-slab 
geometry with p — 1/4 should yield a good approxima- 
tion to F(x,0) in film geometry near bulk T c . 
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FIG. 1: Amplitude of the scaling function of the free energy 
density F(0,p) at T c for n = oo (Eq. <(3j , thin solid line) 
and n = l (Eq. ©, thick solid and dotted lines) in three 
dimensions as a function of the aspect ratio p = L/L\\. MC 
result (full circle) from Ref. [3] for the d — 3 Ising model in 
a cube, p = 1. 

One expects that a similar situation holds for n = 1. 
This is indeed supported by our analytic prediction for 
F(Q, p) as presented in Eq. © below which is shown in 
Fig. 1 as thick solid and dotted lines. We have derived 
this result on the basis of an improved version [j| of the 
lowest-mode separation approach [3, EE EH • Our result 
for F(0,p) agrees very well with earlier MC data [l8| in 
a cubic geometry at p — 1 (full circle in Fig. 1 ). In the 
range p > 1/4 (thick sold line), F(0, p) has the expected 
negative slope. In the range p < 1/4, however, the posi- 
tive slope of the dotted portion of the curve constitutes 
a clear indication for the expected deterioration of the 
quality of the lowest-mode separation approach. In such 
a flat geometry with L/L\\ < 1/4 the system is already 
close to film geometry such that the higher modes are not 



well separated from the single lowest mode. On the other 
hand, together with the result for n = oo, Fig. 1 suggests 
that a calculation of F(x, p) for n — 1 at p — 1/4 should 
yield an acceptable approximation to F(x,0) in film ge- 
ometry near bulk T c . 

Our derivation of F(x, p) for n = 1 is based on the 
Hamiltonian |1| where the decomposition tp = <I> + a is 
made into a homogeneous lowest-mode amplitude $ and 
higher-mode fluctuations a. After integration over a, the 
free energy density is obtained in the form 



f = fo-V- 1 hx / d$exp[-iJ (*)-r(*)] (5) 



with H a ($) = V(\r § 2 + u a ^ 4 ) where f is indepen- 
dent of ro and uq. The higher- mode contribution T($) 
is calculated in one-loop order and expanded around the 
lowest-mode average M 2 = J d$<I> 2 e _ffo / J d$e~ H ° up 
to 0(($ 2 — Mq) 2 ). In truncating our expansion of r($) 
we require that, in the central finite-size region includ- 
ing T — T c , terms of 0(v,q ) are neglected. As we 
are working at fixed dimension 2 < d < 4 there is no 
necessity of further expanding the exponential function 
e~ Ha - T . Thus we maintain the exponential structure of 
the integrand in ([5|) . The resulting bare perturbation ex- 
pression for / contains the bare bulk free energy density 
= limy^oo / in one-loop order above (+) and below 
(— ) T c . The dependence of / — on the aspect ratio p 
appears (i) on the level of the lowest-mode Hamiltonian 
Ho(Q) and (ii) on the level of the contribution of T($). 
The former dependence (i) comes from M§(ro, L, p) — 
(LV~ d wor 1/2 ?%o) with y = ro^V^/wo) 172 where 



■&(y) = / dzz 2 e-i yz - z / \ dze-* yz 
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The latter dependence (ii) is contained in the difference 
between sums over higher modes and bulk integrals in 
wave vector (k) space such as 
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The bare perturbation result needs, of course, to be 
renormalized. Within the minimal renormalization 
scheme in three dimensions 11311 we have obtained the 
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following scaling function for n = 1 
I 3 
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where 



F(x,p) = l^p-^ATTu*)- 1 / 2 ^2Au*a(x, p)R 2 (l, p) 
+Q*xr 1/v l + X8u*R 2 (l,p) + 12u*R 1 (l,p)\ (10) 

with a(x,p) = 12u* 1/2 /- 3 / 2 j o7r 1 / 2 ,9( 2/ ) and 6(ai,p) = 
3? 1 / 2 M* 1/2 / 97r- 3 / 2 z9(y), w* = 0.0412, Q* = 0.945, v = 
(2 — a)/3 = 0.6335. The functions l(x,p), y(x,p), and 
Ri(l,p) are determined by 

y + 12?%) = (0 - 1 / 3 / 2 (47r U *)- 1 / 2 , (11) 
2/ = x Q* r a/(2,/) p- 1 (47Tu*)- 1/2 , (12) 
p) = 4tt(1 - p 2 )r 3 + (/^)- 1 / 1 (/ 2 ,p), (13) 



fl 2 (7,p) = -- + 4tt(1 - p 2 )/- 3 + (Z/4tt 3 )J 2 (/ 2 , p) . (14) 

For hnite p > 0, F(x, p) is an analytic function of x near 
x = 0, in agreement with general analyticity require- 
ments. 

The singular part of the bulk free energy density 
ff b {t) = A ± \t\ d " is, of course, independent of p. It can 
be written as ff b = L~ d F b (x) where F b (x) is the bulk 
part of F(x, p) which is obtained from @ in the limit of 
large \x\. It is given by 



(A-/A+)Q 1 | x 



\du 



for T>T C1 
for T < T c , 



(15) 



with universal numbers Q\ = —0.119 and A~ /A + = 2.04 
in three dimensions. Thus the singular part of the excess 
free energy density f^ x = f s — ff b has the scaling form 



f!*(t,L,L l{ ) = L- d F™(x,p), 
F^(x,p)=F(x,p)-F b ± (x). 



(16) 
(17) 
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FIG. 2: Scaling function F ex (x,p) of the excess free energy 
density for n = l, Eqs. (EJ) , (JTlIl , (U7J) for p = 1/2, 1/3, 1/4 in 
three dimensions as a function of the scaling variable x. 



The prediction of the p dependence of F ex (x,p) above, 
at, and below T c for the Ising universality class as de- 
scribed by (|9|)- (fT7|) without any adjustment of param- 
eters is the central result of this paper. This function 
contains a p dependent minimum slightly below T c . The 
scaling function F ex (x,p) is shown in Fig. 2 for several 
values of p. As expected on the basis of Fig. 1, the differ- 
ence between F ex (x, p) for p = 1/3 and p — 1/4 is rather 
small. 
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FIG. 3: Scaling function X(x,p), Eqs. Q, {H), (QD of 

the Casimir force in three dimensions for finite-slab geometry: 
tp 4 theory at d — 3 (this work, solid line) for p — 1/4, MC 
data [3 (stars) for p = 1/6. 

We note that so far no confirmation of the theory at 
p = 1 j| by MC simulations has been presented except 
right at T = T c [l8| . In particular the prediction of a 
minimum of F(x, 1) below T c for p = 1 [5j is as yet un- 
confirmed since no MC data are available as of yet in 
this regime. For this reason it is particularly interest- 
ing to present here our prediction of a minimum of the 
Casimir force scaling function for small p slightly below 
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T c which can be compared with recent MC data 0, [l(| ■ 
We define the critical Casimir force Fc as imir per unit 
area in a finite-slab geometry as 

F Ca H m ir(t,L,L\\) = = L~ d X(x, p) (18) 

where the derivative is taken at fixed Lu. This defini- 
tion is equivalent to its lattice counterpart introduced in 
(To| . The Casimir force scaling function X can then be 
expressed in terms of F ex as 

v ox op 

(19) 
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FIG. 4: Magnified plot of Fig. 3 near T c : MC estimate for 
Xc 1C at T c (triangle) [13], MC estimate for the minimum 
Xmin (circle) [hj, two- and three- loop e expansion results at 
T c (square and diamond) for p = [all]- 



Numerical evaluation of (fl9|) yields the curve shown 
in Fig. 3. There is good overall agreement with the 
MC data of pjj (with p — 1/6 and L — 20) and with 
the MC data of Q (not shown in Fig. 3) in the range 
~ 2 < x < 15. Somewhat unexpectedly, our result ex- 
hibits a small shoulder near T c . This shoulder is not 
present in the scaling function of the excess free energy 
density F ex (x,p), Fig. 2, but arises through the deriva- 
tive term — (x/v)dF ex (x,p)/dx. 

A more detailed comparison with earlier results is 
shown in Fig. 4. Most significant is the satisfactory 
agreement of the position of the minimum of the theo- 
retical curve x m i n — —0.715 with the MC estimate [Tol | 
x min — —0.681 (full circle in Fig. 4). There is also reason- 
able agreement with regard to the depth of the theoret- 
ical minimum X(x m i ni 1/4) = —0.350 compared to the 
MC estimate X™g = -0.329 (full circle in Fig.4). 
Furthermore, our result X(0, 1/4) = -0.326 at T c is in 
substantially improved agreement with the MC estimate 
[l(| Xf c = -0.304 at T c (triangle in Fig. 4), compared 
to the earlier e expansion results —0.221 in two-loop or- 
der and —0.393 in three-loop order [§] (shown in Fig. 
4 as square and diamond, respectively). 

In summary, we have presented a new approach to the 
analytic calculation of the critical Casimir force scaling 
function in slab geometry for isotropic systems in the 
Ising universality class and have obtained quantitative 
agreement with MC data for periodic boundary condi- 
tions. This approach can be extended to realistic bound- 
ary conditions and to other universality classes which 
may then lead to a satisfactory explanation of the mini- 
mum of the critical Casimir force scaling function below 
T c in real systems [111 ]. 
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